Constraints on the Bulk Lorentz Factor of GRB 990123 
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Abstract. GRB 990123 was a long, complex gamma-ray burst accompanied by an extremely bright optical 
flash. We present the collective constraints on the bulk Lorentz factor for this burst based on estimates from 
burst kinematics, synchrotron spectral decay, prompt radio flash observations, and prompt emission pulse width. 
Combination of these constraints leads to an average bulk Lorentz factor for GRB 990123 of To = 1000 ±100 
which implies a baryon loading of M 
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X 10 Mp. We find these constraints to be consistent with 



the speculation that the optical light is emission from the reverse shock component of the external shock. In 
addition, we find the implied value of Mjet to be in accordance with theoretical estimates: the baryonic loading is 
sufficiently small to allow acceleration of the outflow to F > 100. 
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INTRODUCTION 

The discovery of a prompt and extremely bright opti- 
cal flash in GRB 990123 (Akerlof et al. 1999), imply- 
ing an apparent peak (isotropic) optical luminosity of 
5 X lO'*^ erg s^', has lead to widespread speculation that 
the observed radiation arose from the reverse shock com- 
ponent of the burst. The reverse shock propagates into 
the adiabatically cooled particles of the coasting ejecta, 
decelerating the shell particles and shocking the shell 
material with an amount of internal energy comparable 
to that of the material shocked by the forward shock. 
The typical temperature in the reverse-shocked fluid is, 
however, considerably lower than that of the forward- 
shocked fluid. Consequently, the typical frequency of 
the synchrotron emission from the reverse shock peaks 
at lower energy. It is believed to account for the bright 
prompt optical emission from GRB 990123. The re- 
verse shock emission stops once the entire shell has been 
shocked and the reverse shock reaches the inner edge of 
the fluid. Unlike the continuous forward shock, the hy- 
drodynamic evolution of the reverse-shocked ejecta is 
more fragile. As we will demonstrate, the temperature 
of the reverse-shocked fluid is expected to be mildly- 
relativistic for GRB 990123 and thus the evolution of the 
ejecta deviates from the Blandford & McKee (BM) so- 
lution that determines the late profile of the decelerating 
shell and the external medium. In general, the determi- 
nation of Fq for optical flashes, would play an important 
role in discriminating between cold and hot shell evolu- 
tion. Moreover, the strong dependence of the peak time 
of the optical flash on the bulk Lorentz factor F provides 
a way to estimate this elusive parameter We discuss how 



different constraints on Fq for GRB 990123 are consis- 
tent with optical emission from the reverse shock. We as- 
sume Hq = 65 kms^' Mpc^', Q.yi — 0.3, and Q./^ = 0.7. 



THE ROLE OF T 

Relativistic source expansion plays a crucial role in vir- 
tually all current GRB models (Piran 1999; Meszaros 
2001). The Lorentz factor is not, however, well deter- 
mined by observations. The lack of apparent photon- 
photon attenuation up to « 0.1 GeV implies only a lower 
limit F « 30 (Meszaros, Laguna & Rees 1993), while the 
observed pulse width evolution in the gamma-ray phase 
eliminates scenarios in which F >> 10^ (Lazzati, Ghis- 
elHni & Celotti 2001; Ramirez-Ruiz & Fenimore 2000). 
The initial Lorentz factor is set by the baryon loading, 
that is ifiQC^, where mo is the mass of the expanding 
ejecta. This energy must be converted to radiation in 
an optically-thin region, as the observed bursts are non- 
thermal. The radius of transparency of the ejecta is 
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where E is the isotropic equivalent energy generated by 
the central site. The highly variable y-ray light curves can 
be understood in terms of internal shocks produced by 
velocity variations within the relativistic outflow (Rees 
& Meszaros 1994). In an unsteady outflow, if F were to 
vary by a factor of « 2 on a timescale 5r, then internal 
shocks would develop at a distance Ri w F^c5r > R^^. 
This is followed by the development of a blast wave 



expanding into the external medium, and a reverse shock 
moving back into the ejecta. The inertia of the swept-up 
external matter decelerates the shell ejecta significantly 
by the time it reaches the deceleration radius (Meszaros 
& Rees 1993), 

/ E \ '/^ 
Ra=[ 5^ ^lO'^-lO^cm. (2) 

Given a certain external baryon density no, the initial 
Lorentz factor then strongly determines where both in- 
ternal and external shocks develop. Changes in Fq will 
modify the dynamics of the shock deceleration and the 
manifestations of the afterglow emission. 

REVERSE SHOCK EMISSION AND Tq 

It has been predicted that the reverse shock produces a 
prompt optical flash brighter than 15th magnitude with 
reasonable energy requirements of no more than a few 
10^^ erg emitted isotropically (Meszaros & Rees 1999; 
Sari & Piran 1999a; hereinafter MR99 and SP99). The 
forward shock emission is continuous, but the reverse 
shock terminates once the shock has crossed the shell and 
the cooling frequency has dropped below the observed 
range. The reverse shock contains, at the time it crosses 
the shell, an amount of energy comparable to that in the 
forward one. However, its effective temperature is signif- 
icantly lower (typically by a factor of F). Using the shock 
jump conditions and assuming that the electrons and the 
magnetic field acquire a fraction of the equipartition en- 
ergy Ze and Zb respectively, one can describe the hydro- 
dynamic and magnetic conditions behind the shock. 

The reverse shock synchrotron spectrum is determined 
by the ordering of three break frequencies, the self- 
absorption frequency Va, the cooling frequency and 
the characteristic synchrotron frequency \,„, which are 
easily calculated by comparing them to those of the for- 
ward shock (MR99; SP99; Panaitescu & Kumar 2000). 
The equality of energy density across the contact discon- 
tinuity suggests that the magnetic fields in both regions 
are of comparable strength. 

Assuming that the forward and reverse shocks both 
move with a similar Lorentz factor, the reverse shock 
synchrotron frequency is given by 

v,n = 5.840 X 10'3e2_ie^/!2«aorg,2(l +z)"'/^Hz, (3) 

while the cooling frequency is equal to that of the 
forward shock. Here we adopt the convention Q — IQ^ Qx 
for expressing the physical parameters, using cgs units. 
The spectral power /v,„ at the characteristic synchrotron 
frequency is 

Fv„ = 4.17D2tel(!2^53niJro^2(l +zfhy. (4) 



The distribution of the injected electrons is assumed to 
be a power law of index —p, above a minimum Lorentz 
factor Y,. For an adiabatic blast wave, the correspond- 
ing spectral flux at a given frequency above v,„ is ~ 
Fv„{y /Vm)^^P^^^^^, while below v,„ is characterised by 
a synchrotron tail with Fy ~ ^v,„ (v/v,n)'''^. Similar re- 
lations to those found for a radiative forward shock hold 
for the reverse shock (Kobayashi 2000; hereinafter KOO). 

Unlike the synchrotron spectrum, the afterglow light 
curve at a fixed frequency strongly depends on the hydro- 
dynamics of the relativistic shell, which determines the 
temporal evolution of the break frequencies v,„ and v^. 
The forward shock is always highly relativistic and thus 
is successfully described using the relativistic generali- 
sation of theory of supernova remnants. In contrast, the 
reverse shock can be mildly relativistic. In this regime, 
the shocked shell is unable to heat the ejecta to suffi- 
ciently high temperatures and its evolution deviates from 
the BM solution (Kobayashi & Sari 2000; hereinafter 
KSOO). Shells satisfying 

^^OmE',i%;%i%f>l, (5) 

are likely to have a Newtonian reverse shock[], otherwise 
the reverse shock is relativistic and it considerably de- 
celerates the ejecta. The width of the shell. A, can be 
inferred directly from the observed burst duration by 
A = c7dui/(l +z) assuming the shell does not undergo 
significant spreading (Piran 1999). 

If ^ > 1 , then the reverse shock is in the sub-relativistic 
temperature regime for which there are no known ana- 
lytical solutions. In order to constrain the evolution of F 
in this regime it is common to assume F oc R^s where 
3/2 < ^ < 7/2 (MR99; KSOO). For an adiabatic expan- 
sion. Foe 7-g/(i+2g) and so v,„ oc 7^-3(8+5g)/7(i+2g) and 
Fv„, °= r-('2+iig)/7(i+2g)^ rpj^g spectral flux at a given fre- 
quency expected from the reverse shock gas drops then as 
7^-2(2+3jr)/7(i+2g) i^giQ^ and r-(^+24p+i5p^)/i4(i+2;r) 

above. For a typical spectral index p = 2.5, the flux decay 
index varies in a relatively narrow range (« 0.4) between 
limiting values of g. 



CONSTRAINTS ON THE To OF GRB 
990123 

Despite ongoing observational attempts, the optical flash 
associated with GRB 990123 remains the only event of 
its kind detected to date. Observations of this optical 



It should be remarked that equations (3) and (4) refer to this reverse 
shock regime. Equations for the relativistic case are relatively similar, 
the biggest discrepancy being that the peak flux is inversely propor- 
tional to r (see equations (7)-(9) of KOO). 




FIGURE 1. BATSE and ROTSE light curves for GRB 
990123 as a function of time from the BATSE trigger. Dashed 
lines represent theoretical predictions for the rise r3p-3/2 
and decay oc j—(2i+7ip)/96 adiabatic reverse shock light 
curve, assuming the shell is thin and cooling slowly. We predict 
that the optical flash peaked 41 ± 6 seconds after the trigger. 



flash appear to be in good agreement with early pre- 
dictions for the reverse shock emission (Sari & Piran 
1999b). As a result, numerous studies have been done on 
this event in which reverse shock theory has been applied 
to burst observations in order to constrain physical pa- 
rameters and burst properties, including Fq. Current esti- 
mates on the bulk Lorentz factor for GRB 990123 stretch 
over nearly an order of magnitude, with values ranging 
from « 200 (SP99) to « 1200 (Wang, Dai & Lu 2000). 
It should be noted, however, that these estimates were 
made before accurate burst parameters for GRB 990123 
were known, and consequently they include approxima- 
tions and parameters from other GRB afterglows. 

By fitting multi-frequency afterglow light curves, 
physical parameters for 8 GRBs have recently been re- 
ported (Panaitescu & Kumar 2001b, hereinafter PKOl). 
Best fit values presented for GRB 990123 are 50 
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a rough estimate for the bulk Lorentz factor of Fq = 
1400 ±700 (Panaitescu & Kumar 2001a). Using these 
physical parameters, we present a comprehensive exami- 
nation of the constraints on Fq and report a best-fit value 
based on an analysis of these constraints. 

Time of Peak Flux. Observational estimates for the 
time of peak flux enabled a measurement of the initial 
Lorentz factor with reasonable accuracy using the phys- 
ical parameters specific to GRB 990123. Assuming the 
optical flash was the result of the reverse shock, the ini- 
tial bulk Lorentz factor 
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where 7d is the time of peak flux in the observer frame. 
Light curves for the optical flash and y-ray emission 
are shown in Figure 1. Dashed lines represent theo- 
retical predictions for the rise «: 7^3^-3/2 j^jj^j jgcay 
y--(2i+73p)/96 adiabatic reverse shock light curve, 
assuming the shell is thin, i.e. A < [E / (Inom^c^TQ))^!^ , 
and cooling slowly (KOO). Using the recently reported 
physical parameters, one finds GRB 990123 to have a 
marginal thickness, as predicted by KOO. The observed 
rise time is, however, in good agreement with that of 
a thin shell ~ T^-^ (in contrast with « r'/^ for a thick 
shell). The shape of the light curve is determined by the 
time evolution of the three spectral break frequencies, 
which in turn depend on the hydrodynamical evolution 
of flie fireball. In flie case of GRB 990123, the typical 
synchrotron frequency v,„ = 1.5 x lO'^ Hz is well be- 
low the cooling frequency Vc = 1.0 x lO'^ Hz, and there- 
fore places the burst in a regime with a flux decay gov- 
erned by the relation Ej <^ 7'-(2i+73p)/96 jj^jg implies a 
decay of « T'^ for the optical flash of GRB 990123. 
Applying these light curve predictions to the prompt op- 
tical data and taking observational uncertainties as well 
as burst parameter uncertainties into account, we predict 
that the optical flash peaked Tpoak ~ 41 ±6 s after the 
GRB started. Substitution for Jpeak in equation (6) gives 
Fo = 770±50. 

Synchrotron Spectral Decay. Observations of the op- 
tical peak brightness enable further accurate constraints 
on the value of Fo. The synchrotron spectrum from rel- 
ativistic electrons comprises four power-law segments, 
separated by three critical frequencies. The prompt opti- 
cal flash in GRB 990123 is observed at a frequency that 
falls well below the coofing frequency, but above the typ- 
ical synchrotron frequency: Vq < < Vobs < ^c- The 
synchrotron spectrum for this spectral segment is given 
by F,h. = ^'vJv„,,,/v,„)(''-')/2(i +^)i/2-p/8 ^here v^,, 

is taken to be the ROTSE optical frequency. Assuming 
the optical peak flux Fobs observed in GRB 990123 is 
radiation arising from the reverse shock, we find Fq = 

Radio Flare Observations. In addition to the optical 
flash associated with GRB 991023, a short (-- 30 hr.) ra- 
dio flare was also observed (Kulkarni et al., 1999). Since 
the emission did not grow stronger with time thereby 
demonstrating the properties of a typical radio afterglow, 
it was subsequently proposed that both the prompt op- 
tical emission and the radio flare arose from the reverse 
shock (SP99). As the reverse shock cools the emission 
shifts to lower frequencies which rapidly approach the 
observational radio band. The received flux peaks when 
Vm crosses the band while decreasing mildly relativis- 
tically as v„ oc r-3(8+5g)/7(i+2g) At ? < 1.2 days after 
the initial GRB trigger, the radio flare from 990123 was 
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FIGURE 2. Collective constraints on To for GRB 990123. 
These include estimates from the burst kinematics (narrowest 
distribution), synchrotron spectral decay (widest distribution), 
radio flash observations (medium width distribution), prompt 
emission pulse width (filled arrow), and jet modelling (unfilled 
arrow). We find a best fit value of To = 1000 ± 100. 



observed to be increasing thereby implying Vm > Vobs- 
Using eqn (3) and the evolution of Vm, we compute the 
To which places the peak frequency in the radio band at 
f ~ 1 .2 days. We find Tq = 2000+^{J5^. 

Prompt Emission Pulse Width. Although observa- 
tions of a reverse shock induced optical peak enable 
fairly accurate calculations of the bulk Lorentz factor, 
it remains possible to obtain information on To in situ- 
ations when an optical flash has not been detected. Con- 
sider an internal shock which produces an instantaneous 
burst of isotropic y-ray emission at a time, f, and ra- 
dius, Ri, in the frame of the central engine. The kine- 
matics of colliding shells implies that although photons 
are emitted simultaneously, the curvature of the emit- 
ting shell spreads the aiTival time of the emission over 
a period of ATp, thereby producing the observed width 
in individual pulses. The delay in aiTival time between 
on-axis photons and those at 9 « Ms a function of 
the radius of emission and the Lorentz factor according 
to: Arp/(1 +z) = Ri/{2cr^) where T = To in the early 
phase of the expansion (Fenimore, Madras & Nayashin 
1996). In order to allow photons to escape, Ri must 
be larger than the radius of transparency Rx- This im- 
poses a lower limit on the initial bulk Lorentz factor 
such that: To > (/?x/2cArp)'/2(l +z)'/^. We determine 
ATp for GRB 990123 by measuring the average pulse 
width through autocoiTelation methods based on those 
described in Fenimore, Ramirez-Ruiz & Wu (1999) and 
find ATp « 0.45 s. Using the burst parameters to estimate 
Rx from equation ( 1 ) and applying the inequality relation 
defined above, we find a lower limit of Fq > 200. Fig- 
ure 2 displays the collective constraints on Fq for GRB 
990123. An additional lower limit of Fq > 300 constraint 



from afterglow modelling is also included (Panaitescu & 
Kumar 2002). The combination of these constraints leads 
to an average bulk Lorentz factor for GRB 990123 of 
Fo = 1000 ± 100 which impUes 7?^ « 1.3 x lO''* and a 
baryon loading of Mjei = 8+2^ x lO^^M©. 



CONCLUSIONS 

We have shown that the collective constraints for the 
bulk Lorentz factor of GRB 991023 are compatible with 
current reverse shock theory. In addition, our best fit 
value of Fq ~ 1000 provides confirmation of the ultra- 
relativistic nature of GRBs. The implied values for R^ 
and Mjer are in accordance with GRB theory: the radius 
of transparency is within theoretical estimates and the 
baryonic loading is sufficiently small to allow acceler- 
ation of the outflow to F > 100. As we have discussed, 
the bulk Lorentz factor plays a crucial role throughout 
all stages of GRB evolution, and therefore the ability to 
constrain Fq provides clues on the nature of gamma-ray 
bursts. 
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